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1.1 Introduction 

Environmental and genetic mutations can transform the cells in a co-operating 
healthy tissue into an ecosystem of individualistic tumour cells that compete 
for space and resources [TJ |2J [3] . If we consider a tumour as an ecosystem it 
is possible to utilise tools traditionally used by ecologists to study the evolu- 
tion of a population in which there is some degree of phenotypical diversity. 
One such tool is evolutionary game theory (EGT) which merges traditional 
game theory with population biology gj . It allows the prediction of successful 
phcnotypes and their adaptation to environmental selection forces. EGT is 
considered as a promising tool in which to frame oncological problems [5] and 
has been recently made more relevant by phenotypic studies of carcinogenesis 
such as the ones by Hanahan, Weinberg and colleagues [ElE]. 

Game theory (GT) was introduced by von Neumann and Morgenstern as 
an instrument to study human behaviour [5J [3] • A game describes the inter- 
actions of two or more players that follow two or more well defined strategies 
in which the benefit of each player (payoff) results from these interactions 
PU| . GT can be employed to study situations in which several players make 
decisions in order to maximise their own benefit. GT was initially introduced 
to model problems in economics, social and behavioural sciences and is used 
as a formal way to analyse interactions between agents that behave strategi- 
cally. Evolutionary game theory is the application of conventional GT as used 
by economists and sociologists to study evolution and population ecology [4]. 
As opposed to conventional GT, in EGT the behaviour of the players is not 
assumed to be based on rational payoff maximisation but it is thought to 
have been shaped by trial and error - adaptation through natural selection or 
individual learning [11 . In the context of the evolution of populations there 
are two GT concepts that have to be interpreted in a different light. First, 
a strategy is not a deliberate course of action but a phenotypic trait. The 
payoff is Darwinian fitness, that is, average reproductive success. Secondly, 
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the players are members of a population that compete or cooperate to obtain 
a larger share of the population 0]. 

To illustrate some of the ideas in EGT let us consider the following example 
named the Hawk-Dove game In this game we study an imaginary popu- 
lation of individuals and a resource V which affects the reproductive success 
of the individuals in this population. The population contains two phenotypes 
that represent two different strategies to access the resource. When two indi- 
viduals compete for the resource the outcome will depend on the phenotypic 
strategies involved. The first phenotype, called Hawk in the game, always es- 
calates the fight until injured (at a cost in fitness equal to C) or until the rival 
retreats. The second phenotype, known as Dove in the game, will retreat if the 
opponent escalates, that is, if the opponent seems determined to fight. The 
interactions between the different phenotypes are shown in the payoff table 

o 

Table 1.1. Payoff table for the change in fitness in the Hawk-Dove game. 





Hawk 


Dove 


Hawk 


v-c 
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Dove 


V 


V 
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Table |1.1| presents the interactions between the different phenotypes con- 
sidered in the game. The table should be read following the columns such that 
the payoff for a Hawk playing another Hawk is ~ expressing the fact that 
they both have to share the resource and that they stand an equal chance 
of getting injured. The payoff of a Hawk playing a Dove is V since the Dove 
will withdraw from the competition. A Dove playing a Hawk gets no payoff 
since it withdraws and when playing another Dove it will get the resource V 
in half of the occasions. With this information it is possible to predict that 
if the population is mainly composed of individuals with the Dove phenotype 
then a Hawk individual will have a significant fitness advantage (as in most 
of the interactions, the rivals are likely to be a Dove and thus retreat from 
a full scale fight for the resource). On the other hand if the fitness cost of 
injury is more than twice as high as the benefit provided by the contested 
resource then a Dove would be quite successful in a population dominated by 
Hawks (since a Hawk that interacts frequently with other Hawks is likely to 
be eventually wounded and a Dove will always avoid costly wounds). Using 
this example Maynard Smith introduced the concept of an evolutionary stable 
strategy (ESS) [IT]. An ESS is defined as a phenotype that, if adopted by the 
vast majority of a population, will not be displaced by any other phenotype 
that could appear in the population as a result of evolution [TT] . Under this 
definition the Dove phenotype cannot be an ESS and only under some specific 
circumstances (when the fitness benefit of getting the resource outweighs the 
fitness cost of an injury) would a Hawk phenotype be evolutionary stable. 
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GT has been used to address many problems in biology in which different 
species or phenotypes within one species compete. Examples of this are the 
evolution of sex ratios [H] , the emergence of animal communication [T3] and 
fighting behaviour and territoriality [11J. A recent focus on the capabilities 
that cells acquire as tumours evolve [5] has shown how the interplay between 
different phenotypes with different capabilities can lead to different evolution- 
ary paths. This stresses the importance of GT as a modelling tool in cancer. 
The most important capabilities which cells have to acquire in a neoplasm 
that will become a malignant tumour are shown in figure |1.1| They include: 
unlimited replicative potential, environmental independence for growth, eva- 
sion of apoptosis, angiogenesis and invasion. The circumstances in which these 
capabilities evolve and spread through the tumour population can be studied 
using GT. Some of the most important milestones in the transformation of a 
healthy tissue cells into malignant cancer such as tumourigenesis, angiogenesis 
and invasion, have already been approached with GT. 
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Fig. 1.1. Acquired capabilities of cancer: Hanahan and Weinberg suggest that most 
if not all cancers have acquired the same set of capabilities during their evolution. 
These capabilities are: unlimited replicative potential, environmental independence 
from anti growth signals, production of own growth signals, evasion of apoptosis, 
angiogenesis and invasion (from [S] with permission). 



The remaining of this chapter will provide, to the best of our knowledge, 
all the relevant examples of the application of GT to the study of the somatic 
cancer evolution and finally hint some of the possible future venues of this 
method in the context of cancer research. 
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1.2 Tumourigenesis 

Tumour initiation requires the acquisition of a number of phenotypic capabili- 
ties such as evasion of apoptosis and independence from environmental signals 
(see figure The evolution of these capabilities, normally acquired when 
the tumour is still in the avascular stage, are studied in research by Tomlinson 
and Bodmer OfTS] and by Gatenby and Vincent [15] . 

1.2.1 Evasion of apoptosis 

Problem. Apoptosis or programmed cell death is a mechanism that hinders 
tumour progression. Cells with a working apoptotic machinery die when ge- 
netic abnormalities are detected [6 . Thus cells in a malignant cancer have to 
evolve mechanisms to disable the apoptotic machinery. 

Model. Tomlinson and Bodmer [TJ] present a model in which three dif- 
ferent apoptosis evasion related strategies are considered: 

1. Cells that produce a paracrine growth factor to prevent apoptosis of neigh- 
bouring cells. 

2. Cells that produce an autocrine growth factor to prevent apoptosis of 
themselves. 

3. Cells susceptible to paracrine growth factors but incapable of production 
of factors. 

The aim of the model is to study the possibility of stable coexistence of 
the different phenotypes (polymorphism) that could be possible in a tumour 
when only these three phenotypes are considered. 



Table 1.2. Payoff table for the game of programmed cell death. Parameter a rep- 
resents the cost of producing a paracrine factor, b is the cost of the factor produced 
in an autocrine fashion, c is the fitness benefit of evading apoptosis. 





1 


2 


3 


1 


1 - a + b 


l + b + c 


1 + 6 


2 


1 - a 


l + c 


1 


3 


1 - a 


1 + c 


1 



In table 1.2 a is the cost of producing the paracrine factor, b the benefit 
of receiving the paracrine factor and c the benefit of producing the autocrine 
factor. 

Results. If a is positive then the third strategy displaces the first one 
from the population. If only the other two strategics arc considered then if the 
benefit of the autocrine factor, provided by c, is positive the second strategy 
will displace the third one. In most relevant situations the model shows a 
strong selection for the autocrine factor producing phcnotype and under the 
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assumptions of the model the altruistic strategy (the first one) will always be 
displaced. 

Remarks. The model is very simple and easy to understand while at the 
same time it captures the relevant features necessary to study the evolution of 
the mechanisms to avoid apoptosis. However the authors do not explain the 
mechanisms by which these phenotypes could appear in a tumour and what 
would be the biological explanation of the costs and benefits of the different 
growth factors. 

1.2.2 Environmental poisoning 

Problem. Tomlison introduced a further model in which he considers the 
hypothesis that tumour cells might boost their own replicative potential at 
the expense of other tumour cells by evolving the capability of producing 
cytotoxic substances [T5] , 

Model. Tomlinson speculates with different strategies that cells may 
adopt to produce or cope with toxic factors. The main model aims to study 
the polymorphic equilibria when cells can adopt one of the three following 
strategies: 

1. Cells producing cytotoxic substances against other cells, 

2. cells producing resistance to external cytotoxic substances, and 

3. cells producing neither cytotoxins nor resistance. 

Table [T73| shows the payoff table of the game with these phenotypes. 

Table 1.3. Payoff table for the change in fitness for cells in a tumour in which the 
base payoff is z, e the cost of producing the cytotoxin, f the fitness cost of being 
affected by the cytotoxin, g the advantage of subjecting another cell to the cytotoin 
and finally the cost of developing resistance to the cytotoxin is h. 
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Results. Game theoretical analysis and simulations show that production 
of cytotoxic substances against other tumour cells can evolve in a tumour 
population and that several cytotoxin related strategies may be present at a 
given time (polymorphism). 

Remarks. Although the author admits that there is little experimental 
evidence for mutations that cause tumour cells to harm their neighbours, the 
Warburg effect could fit nicely in the framework presented in this work. The 
Warburg effect describes the switch of tumour cells from the conventional 
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aerobic metabolism to the glycolytic metabolism. This metabolism is less ef- 
ficient but produces, as a by-product, acid that can harm neighbouring cells 
[T71 [T5J Q20 [20] ■ Thus, in the game described by Tomlinson, e could correspond 
to the fitness loss of the less efficient glycolytic metabolism, f could be the fit- 
ness loss of a normal cell in an acid environment and g the fitness benefit 
received by glycolytic cells that can take advantage of the harm done to their 
non glycolytic neighbours. 

1.3 Angiogenesis 

Problem. A very important capability that has to be acquired by tumours on 
the path to cancer is angiogenesis. Without access to the circulatory system 
tumours do not grow to sizes bigger than 2mm in diameter [5T] . Cells capable 
of angiogenesis produce growth factors that promote the creation of new blood 
vessels that can provide nutrients and oxygen to previously unreachable areas 
in a growing tumour. Presumably the factors will be produced at a cost to 
the tumour cell. 

Model I. In their interpretation of an angiogenic game Tomlinson and 
Bodmer [13] consider two strategies: cells denoted as A+ can produce an- 
giogenic factors at a fitness cost i and cells denoted as A- that produce no 
angiogenic factors. In any case cells will get a benefit j when there is an in- 
teraction involving an angiogenic factor producing cell. The payoffs for the 
interactions between these cells are shown in table 11.41 

Table 1.4. Payoff table for the change in fitness for a cell in a tumour with cells 
capable of producing angiogenic factors (A+) and cells susceptible to benefit from 
growth factors (A-). 





A+ 


A- 


A+ 


1-i+j 


1+j 


A- 


1-i+j 


1 



Results. The model shows that as long as the benefit j of angiogenesis 
is greater than the cost i of producing angiogenic factors then both types of 
strategics will be present in a tumour in proportion to these costs. 

Remarks. This model of angiogenesis is rather simplistic but constitutes a 
nice foundation for later models that take into account spatial considerations. 
A more significant drawback is that the model does not attempt to suggest 
a link between the different fitness costs and benefits and the underlying 
biological mechanisms. 

Model II. The first extension to this model was proposed by Bach et al 
[2"2"| suggesting a game with interactions between three players. In this game 
the benefit j of angiogenesis is obtained only when at least two of the three 
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players produce the angiogenic factor. The new payoff table is shown in table 
IQ1 



Table 1.5. Payoff table for the change in fitness for a cell in a tumour with cells 
capable of producing angiogenic factors (A+) and cells susceptible to benefit from 
growth factors (A-). In this version the interactions involve three players. The payoff 
of a player is given by the columns (ie. the payoff of a A+ cell interacting with a 
A+ and a A- cell is 1 - i + j but a A+ interacting with two A- cells is 1 - j). 





A+ 


A- 


A+, A+ 


1-i+j 


1+j 


A+, A- 


1-i +j 


1 


A-, A- 


1-i 


1 



Results. The authors produce simulations using table 1.5 and show (see 
figure 1.2 1 that this game yields results than differ from the original ones 



introduced by Tomlinson and Bodmer. In this case even when the cost i is 
smaller than the benefit j there are many scenarios for which the angiogenic 
strategy will be displaced from the population. The existence of a polymorphic 
equilibrium containing the angiogenic strategy depends, intuitively, on the 
cost of producing angiogenic factors in comparison to the benefit they give 
but also on the relative frequency of cells playing the angiogenic strategy in the 
population. These results suggest that a gene therapy against the reparation 
of tumour supressor genes would only need to change a fraction of the mutated 
cells before the dynamics of the system drives them to extinction. 

Model III. In a separate research Bach et al [23] studied how a spatial 
version of the angiogenesis game could produce different results to those of 
the original model by Tomlinson and Bodmer. In this game players inhabit a 
100x100 lattice. Each player can follow either the angiogenic (A+) or the non- 
angiogenic (A-) strategy. Time is discrete and in each time step a number of 
cells is removed from the lattice at random. Neighbouring cells compete (using 



table 1.4 1 to occupy unallocated space. The candidate cell that achieves the 



highest score interacting with the neighbours determines the strategy that 



will be followed by the new cell in the vacant slot (see figure 1.3 1. 

Results. The authors found the results of the new formulation of the 
model markedly different from those of the non-spatial counterpart. In spatial 
models, space tends to favour growth promoters in ways that cannot be seen 
in the non-spatial model. In any case polymorphic equilibria do exist and 
in both spatial and non spatial cases the proportion of angiogenic players 
increases as the benefits of angiogenesis increase or as the cost of producing 
the growth factors decreases. The authors also conclude that, contrary to 
other evolutionary models [23] , space does not significantly favour co-operative 
strategies in populations of cells. 



8 David Basanta*, Andreas Deutsch 



a 0.4 



10 15 20 

Generation number 



) 8 6 o o o o a 



Fig. 1.2. Simulations of the threshold model in which the benefit of angiogenesis (j) 
is three times the cost of producing the growth factors (i). The plot shows the results 
of three simulations with a different of A+/A- players in the starting population. 
The ordinate is the frequency of A+ cells and the abscissa is the number of cell 
generations, (from [22] with permission ) 
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Fig. 1.3. A diagrammatic representation of the cell update process in the spatial 
model of the angiogenesis game by Bach and colleagues |23| . First a cell is selected 
to be removed (in this case the C marked in bold in the configuration on the left). 
Then the four neighbouring cells all assess their fitness (middle configuration). The 
cell with the highest fitness then reproduces into the empty space (right), (from [23] 
with permission) 



1.4 Motility /invasion 

A tumour in which cells develop the capability of invading other tissues be- 
comes a malignant tumour and thus significantly worsens the prognosis of a 
patient. EGT is a tool that could greatly improve our understanding of the 
circumstances that influence the successful evolution of invasive phenotypes. 

Model I. Mansury and colleagues have recently introduced a rather un- 
conventional EGT formulation that they used on top of a previous cellular 
automaton model in the context of brain tumours [25] [5B] . This GT module 
allows the original model to deal with cell-cell interactions in a tumour consist- 
ing of cells with different phenotypes. The model was designed to investigate 
the genotype to phenotype link in a polymorphic tumour cell population. 
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This model covers two distinct strategies. Strategy P is characterised by a 
highly proliferative genotype and a high number of gap junctions. Strategy 
M is characterised by a highly migratory phenotype and a low number of 
gap junctions (which are used for cell to cell communication). As opposed to 
other game theoretical models in which a payoff table determines the fitness 
change of players when they interact, in this model three tables are used not 
to compute fitness change but to encode the rules of the cellular automata. 
The tables are shown in tables |1.6| |1.7| and |1.8| 



Table 1.6. Payoff table that describes the rate of change in cell to cell commu- 
nication skills depending on the phenotype of the interacting cell in the model by 
Mansury et al. |26| . 
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Table 1.7. Payoff table that describes how the proliferation capability of a cell 
is influenced by its interaction with cells with other phenotypes (proliferative or 
motile) in the model by Mansury et al. |26j . 
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Table 1.8. Payoff table that describes how the motility of a cell with a given 
phenotype changes depending on what other cells the cell is interacting with. Model 
by Mansury et al. |26| . 
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TT 


M 


TI 
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Table |1.6| shows how the cell to cell communication capabilities change 
according to the strategics followed by the interacting cells. Since the authors 
assume that there is a negative correlation between the extent of communi- 
cation and proliferation then the proliferative capability of a cell with the 
proliferative strategy P will be lower if it interacts with another cell with the 
same strategy. On the other hand it will be higher if it interacts with a cell 
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with the motile strategy M. Table \L7\ shows that cells will have a higher pro- 
liferative potential when they interact with cells that have the motile strategy. 
Table |1.8| shows that interacting with cells with the motile strategy will also 
increase the probability of motility of the interacting cell. These tables are 
used to guide the behaviour of the cells on a 500x500 lattice containing ini- 
tially 5 cells of each type and two unequal sources of nutrients in two different 
locations. In each time step the interactions of a cell with its neighbours are 
used to adjust the probabilities of proliferation and motility. 

Results. The authors used this model to study how varying the payoffs for 
A-A interactions affects a number of features of the tumour such as the speed 
(see figure 1.4 1 at which it reaches the nutrient sources or the fractality of the 
resulting spatial patterns. 




Fig. 1.4. Plot of the time to reach a nutrient source (i.e. the inverse of the tumour's 
average velocity) versus various payoff values conferred to P-P intercellular interac- 
tions (where A is the label for the proliferative phenotype in the paper). The error 
bars indicate the standard deviations from performing 10 Monte Carlo simulations 
for each value of payoffs for P-P interactions, (from [26] with permission). 



Remarks. Although the model is interesting it is not as simple as other 
models reviewed in this chapter. A significant drawback is that the model is 
not evolutionary in the sense that it does not take into account the possibility 
of mutations introducing new phenotypes or tumour cells producing offspring 
different from their progenitors. This has been identified by the authors as 
something to be addressed in a future version of the model although they 
still claim the model to be based on EGT. Moreover, the authors do not take 
advantage of the tools provided by GT analysis to study the different steady 
state situations that could arise for different values of the payoff tables. This 
is probably due to the fact that the model is not a conventional GT model 
and the conventional game theoretical tools would not be easy to use in this 
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context. Also, the model provides a view on the dynamics of the tumour 
growth that is rarely found in most other game theory analysis whose focus 
is on the study of equilibria. 

Model II. The emergence of invasive phenotypes is influenced not only by 
its interaction with one phenotype or the other but by the complex interplay 
of several phenotypes which, in many cases has an indirect effect. Basanta et 
al. |27j hypothesise a number of scenarios in which three types of phenotypes 
interact in games with two and three players. They place these phenotypes in 
the context of an evolutionary non-spatial game theoretical model to test the 
hypothesis by Gatcnby and colleagues that tumour invasion is promoted by 
the emergence of cells with a glycolytic metabolism [20] . The model represents 
a glioma tumour populated by cells with enhanced proliferative capabilities, 
known as autonomous growth cells (AG), which can mutate into cells whose 
phenotype can make then follow either a more motile strategy (INV) or the 
glycolytic metabolism strategy (GLY). 



Table 1.9. Payoff table that represents the change in fitness of a tumour cell with 
a given phenotype interacting with another cell. Three different strategies are con- 
sidered, those with higher replicative potential (AG), enhanced motility (INV) and 
glycolytic metabolism (GLY). The base payoff in a given interaction is equal to 1 
and the cost of moving to another location with respect to the base payoff is c. The 
fitness cost of acidity is n and k is the fitness cost of having a less efficient glycolytic 
metabolism. The table should be read following the columns, thus the fitness change 
for an INV cell interacting with an AG would be 1 — c. 





AG 


INV 


GLY 


AG 
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1 - c 


| + n — k 


INV 
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1 - - 

X 2 


1 - k 


GLY 


\-n 


1 - C 





Results. Table [L9| defines the interactions between the three phenotypes. 
The authors investigate a number of scenarios in which subgames based on 
two strategies are used to study how one phenotype could emerge in a tumour 
populated by cells that use a different strategy. The results show that cells 
with a higher replicative potential (AG) and invasive cells (INV) can coexist 
in a tumour as long as the fitness cost of motility is not too high. They 
also reveal that autonomous growth (AG) cells cannot coexist with glycolytic 
(GLY) cells. An interesting result of the model when all three strategies are 
considered simultaneously is that the appearance of the invasive phenotype 
is facilitated by the existence of glycolytic cells. Figure 1.6 shows how the 
proportion of invasive cells (X axis) increases as the cost of having a glycolytic 
metabolism (k) decreases and the cost of living in an acid environment (n) 
increases. In other words, the success of the invasive phenotype depends on the 
same factors that determine the fate of the glycolytic phenotype. The model 
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suggests that any therapy that could increase the fitness cost of tumour cells 
switching to a glycolytic metabolism or the susceptibility of normal cells to 
acid environments might decrease the probability of the emergence of more 
invasive phenotypes. 




Fig. 1.5. Proportion of invasive cells in a tumour with three phenotypes (au- 
tonomous growth, invasive and glycolytic), k is the cost in terms of fitness of adopting 
the glycolytic metabolism whereas n is the fitness cost of a normal cell when staying 
with a glycolytic cell. 



1.5 Outlook 

The examples described in this chapter treat tumour populations as ecosys- 
tems of potentially co-operating and/or competing cells. In these ecosystems, 
the success of one phenotype depends on its interactions with other existing 
phenotypes. Such an approach has been shown to be a helpful and useful 
way to study cancer evolution. Most of the applications use GT to study the 
steady state of a population of tumour cells that follow different strategies 
dictated by their phenotypes, acquired as a result of genetic and epigenetic 
mutations. This can be very relevant to study the different ways in which 
a tumour population may evolve under different model parameters and as- 
sumptions of their interplay. Such studies could lead to cancer therapies that 
would alter the dynamics of cancer evolution towards benign tumours. One 
limitation of EGT/GT is that all relevant potential phenotypes/stratcgics 
have to be known a priori if a good understanding of cancer evolution is to 
be obtained. Knowing all the potential relevant phenotypes might be difficult 
and even if the phenotypes and their interactions are well known the EGT 
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model might be too complicated to be analysed. Moreover, with more strate- 
gies/phenotypes the composition of the population may not converge to an 
equilibrium and the frequencies of the phenotypes could keep oscillating in a 
regular or chaotic fashion [?]. 

One more limitation of the GT models shown in this chapter is that they 
do not study the dynamics in a tumour population (dynamics which may or 
may not lead to an equilibrium). GT models that make use of population 
biology have the potential to overcome this limitation and also ease the con- 
nection between a quantitative model and experimental data as the payoff 
tables used in more conventional EGT models tend to make the assumption 
that the fitness values are independent of space or time. One promising venue 
is to couple conventional GT with population dynamics which is also based on 
the assumption that successful strategies spread [28, Hj . This trend is shown 
in the work of Gatenby and Vincent [TB] whose EGT model uses methods 
from population biology in order to study how the phenotypes of cells in a 
population evolve towards ESS. Gatenby and Vincent adopted a game the- 
ory approach influenced by population dynamics to study the influence of the 
tumour-host interface in colorectal carcinogenesis. The authors formulated an 
extended system of Lotka-Volterra equations to model the effect of nutrients 
and sensitivity to growth constraints in the proliferation of tumour cells. The 
cells in the tumour population are characterised by the number of substrate 
transporters (a higher number of them allow more nutrients in the cell) on the 
cell surface and by the cell response to normal growth constraints. Initially 
all the tumour cells are assumed to have normal values for both parameters 
but these values are allowed to change and the authors study their evolu- 
tion. Their work demonstrates that normal cells in a multicellular organism 
occupy a ridge-shaped maximum of the fitness landscape that allows the het- 
erogeneous coexistance of multiple cells. This fact makes them susceptible to 
mutations that are fitter and thus allow the somatic evolution that charac- 
terises cancer. The authors also conclude that any therapy that reduces the 
population density would be counterproductive since it would allow the more 
aggressive phenotypes to grow. Subsequent research uses the same formula- 
tion to study other stages of tumour progression like invasion and metastasis 

HOGS]- 

GT is also a suitable tool to frame cooperative effects. Evolution of co- 
operation is a sub- field of GT pioneered by Robert Axelrod [3T] . Researchers 
in this field study the circumstances under which selfish agents will spon- 
taneously co-operate. An example of how evolution of cooperation could be 
used to study cancer evolution was provided by Axelrod and collaborators 
in a recent paper [32J in which they show how tumour cells can cooperate 
sharing skills and capabilities (such as the production of angiogenic factors or 
paracrine growth factors needed to escape the homeostatic regulation of the 
tissue) . 

Evolutionary game theory has a relatively short history in the field of 
theoretical oncology but with an increasingly better understanding of the role 
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Fig. 1.6. (a) Upon reaching equilibrium, normal cells sit on the ridge-shaped maxi- 
mum in the in vivo adaptive landscape. This provides a tissue structure that allows 
coexistence of multiple distinct noncompeting populations, a necessary condition for 
the formation of multicellular organisms, (b) The presence of tumour cells deforms 
the tissue adaptive landscape over time so that normal cells arrive at a local mini- 
mum and are then readly sunject to invasion by the mutant phenotypes (from [16] 
with permission). 

of the microenvironmcnt in tumour evolution [33 and with the recent interest 
in studying cancer from ecological [51 [3] and phenotypic [5] viewpoints, the 
role of this tool to understand the interactions between all the relevant agents 
within a tumour and their role driving cancer evolution will surely rise. 
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